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Abstract 

Using a conjecture that allows to approach separable-variables con- 
ductivity functions, the elements of the Modern Pseudoanalytic Func- 
tion Theory are used, for the first time, to numerically solve the Dirich- 
let boundary value problem of the two-dimensional Electrical Impedance 
Equation, when the conductivity function arises from geometrical figures, 
located within bounded domains. 



1 Introduction 

The study of the solutions corresponding to the Electrical Impedance Equation 

V • (crVw) = 0, (1) 

where a represents the conductivity and u denotes the electric potential, is 
fundamental for the proper understanding of a wide variety of boundary value 
problems, that posses special relevance in different branches of Mathematical 
Physics. Among many important problems, it is particularly interesting the one 
corresponding to the Electrical Impedance Tomography, because of its impor- 
tance in Applied Physics and Engineering. As a matter of fact, only taking into 
account the Medical Imaging applications, no doubts remain about its relevance. 

The study of this inverse problem, correctly posed in mathematical terms 
by A. P. Calderon [3], is based upon iterative methods that employ solutions of 
the forward problem, attempting to fulfil a certain boundary condition (usually 
upcoming from physiological measurements), by introducing variations in the 
conductivity function at every step, in order to reduce the difference between 
the theoretical u and the measured one (see e.g., the classical work [T^). 

Yet, the mathematical complexity of ([T]) posed so strong challenges, that 
the Electrical Impedance Tomography remained only as an alternative Medical 
Imaging technique, seldom considered among the basic clinical applications. As 



a matter of fact, the Electrical Impedance Tomography is still considered an ill 
posed problem. 

In this direction, the mathematical foundations dedicated to this topic, per- 
haps received one of the most important contributions when V. Kravchenko in 
2005 [H], and independently K. Astala and L. Paivarinta in 2006 [T], discovered 
that the two-dimensional case of ([!]) was completely equivalent to a special kind 
of Vekua equation |12j , which had been deeply studied in a variety of works, 
been the most important those published by L. Bers |2j and by I. Vekua [12]. 

The list of novel works arising after such discovering, is as long as it is 
interesting. Still, it shall be pointed out that just some of them are dedicated to 
the Engineering applications, since it is not clear how to model a wide variety 
of physical phenomena, in order to make them susceptible to be analysed by 
means the Modern Theory of Pseudoanalytic Functions (see e.g. [8]). 

This work intends to make a positive contribution in this direction. Already 
in was posed a basic idea for interpolating values of conductivity within 
bounded domains, in order to obtain a certain class of analytic representations: 
Two-dimensional separable-variables functions, one of the keys to fully applied 
the novel mathematical methods into Engineering problems. 

In this work, a more general methodology is posed, based upon the main 
Conjecture arising from Starting with some examples for which the con- 
ductivity functions are known in exact form, the elements of the Pseudoanalytic 
Function Theory are used, for the first time, to approach solutions of the for- 
ward boundary value problem for ([!]), considering conductivity distributions 
upcoming from geometrical figures. 

Even the full set of examples is provided for a circular domain, the results can 
be extended to a wide variety of bounded domains. Thus, the methods provided 
along these pages, could well be ready for studying images corresponding already 
to physical experimental models. 

2 Preliminaries 

Following [2] , let the complex- valued functions F, G satisfy the condition 



where F denotes the complex conjugation of _F: F — ReF — ihaF, and i is the 
standard imaginary unit P = —1. Thus, any complex- valued function W can 
be expressed by means of the linear combination of F and G: 



where 4> and are purely real functions. Two complex functions that fulfil ([2| 
shall be called a generating pair (F, G). Bers [2] introduced the {F, G')-derivative 
of the function W according to the expression: 




(2) 



W ^<I)F + i!G, 



(3) 
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But this derivative will exist if and only if 

(ckcb) F + G = 0. (4) 

Here 

dz=dx-idy, (h^ dx+idy. 

Notice that these operators are classically introduced with the factor i, but it 
will result somehow more convenient to omit it int his work. 
Introducing the functions 

_Fd,G-Gd,F FckG-Gd^F 
FG^GF ' , 
_Fd,G-Gd,F _ Gd^F-Fd^G 

FG-GF ' ^(^^G)-- pG-GF ' 

the expression ([s]) of the G)-derivative will turn into 

d(F.G)W^d,W-A(^F^G)W-B^F^G)W, (6) 

as well the condition Q can be written as 

d^W-a(^F,G)W-b^p^G)W^O. (7) 

The functions defined in ([5]) are called the characteristic coefficients of the 
generating pair {F, G) , and the functions W satisfying the condition ([t]) are 
named (F, G)-pseudoanalytic. Indeed, the equation([7| is know as the Vekua 
equation [12j . and in many senses is the foundation of the present work. 

The following sentences were originally presented in [2] and [8], and they 
have been slightly adapted for the purposes of this paper. 

Theorem 1 The elements oj the generating pair (F, G) are {F, G) -pseudoanalytic: 

d{F.,G)F = d(F,G)G = 0. 

Remark 1 Letp he a non-vanishing function within a hounded domain £7 G M^. 
The functions 

Fo=P, Go = -, (8) 
P 

constitute a generating pair, whose characteristic coefficients are 

^(f'o.Go) = «(Fo,Go) 0' 
^{Fo,Go) = — > HFo,Go) - — • 

Definition 1 Let {Fq, Gq) and {Fi, Gl) be two generating pairs of the form 
and let their characteristic coefficients satisfy the relation 

B(Fi.,Gi) ^ -bi^F^^Go)- (10) 

Thus the pair (Fi,Gi) will be called a successor of the pair {Fq,Go), whereas 
(Fq, Go) will be called a predecessor of [Fi, Gi). 
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Definition 2 Let 

{(F™,G„0}, m = 0,±l,±2,... 

be a set of generating pairs, and let every Gm+i) be a successor of 

{FrmGm)- Thus, the set {{Fm,Gm)} is called a generating sequence. Moreover, 
if there exist a number k such that {Fm,Gm) = {Fm+k,Gm+k) the generating 
sequence is said to be periodic, with period k. 

Finally, if {F,G) = {Fo,Go), the generating pair {F,G) will be embedded 
into the generating sequence {{Fm,Gm)}- 

Theorem 2 Let {F, G) be a generating pair of the form and let p be a 
separable-variables function: 

p = Pi{x)p2{y), 

where G M. Thus {F,G) is embedded into a periodic generating sequence, 
with period k — 2, such that 

P _ P2{y) ^ _ - PiW . 

Pl{X) P2[yj 

when m is even, and 

Frn ^ Pi{x)p2iy), G„ 



Pi{x)P2{.yy 

when m is odd. 

Moreover, if in particular pi{x) = 1, it is easy to see that the generating 
sequence in which {F, G) is embedded will be also periodic, but it will posses 
period k ~ I. 

L. Bers also introduced the concept of the (Fq, Go)-integral of a complex- 
valued function W. The detailed conditions for its existence can be found in 

m- 

Definition 3 Let {Fq,Gq) be a generating pair of the form. Its adjoin gener- 
ating pair {Fq,Gq) is defined according to the formulas 

F* = -iFo, G*Q^-iGo. 

Definition 4 The {Fq,Go) -integral of a complex-valued function W (when it 
exists l^) is defined as: 

j Wd(^Fo,Go)Z ^ FoRe J G*Wdz + GoRe J F*Wdz, 

where T is a rectifiable curve within a domain G C. Specifically, if we consider 
the (Fq^Gq) -integral of the {Fq^Gq)- derivative of W , we will have: 

d(Fo..Go)Wd(^F^^Go)^ = -<i>{zo)F{z) - ^(zo)G(z) + W{z). (11) 

Here z — x + iy, and zq is a fixed point in the complex plane. But according to 
the Theoremfl the {Fq, Gq) -derivative of F and of G vanish identically, thus the 
expression [11] can be considered as the {Fq, GQ)-antiderivative of d(^F^^^Ga)W ■ 
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2.1 Formal Powers 

Definition 5 The formal power (ao, z\ zq) belonging to the generating pair 
(FjmGm), with formal degree (0), complex constant coefficient Oq, center at Zg, 
and depending upon z — x + iy, is defined according to the expression: 

Z^^ (a„, z; zq) = XF{z) + fiG{z), 

where A and fj, are complex constants that fulfil the condition: 

XF{zo) + AiG(zo) oq. 

The formal powers with higher degrees are calculated according to the recursive 
formulas: 

^ (a„, z; zo)=n / Z^'^Zl^ (a„, z; zq) £^(f„,g,„)2, (12) 

J za 

where n = 1,2,3,... Notice the integral operators in the right side of the last 
expression are {Fm,Gm)-antiderivatives. 

Theorem 3 The formal powers posses the following properties: 

1. Every zln' (a„, z; zq) , rt = 0, 1, 2, ... is an {F„i, Gm) -pseudoanalytic func- 
tion. 

2. Let a„ = + ia'^, where and are real constants. The following 
relation holds 

Z(„") (a„, z; zo) - (1, z; z„) + a" Z^^^ {z, z; Zq) . (13) 

Theorem 4 Let W be an {Fm, Gm)-pseudoanalytic function. Then it can be 
expressed by means of the so-called Taylor series in formal powers: 

oo 

W^J2^^nHan,z;zo). (14) 

n=0 

Since any {Fm, Gm) -pseudoanalytic function W accepts this expansion, this is 
an analytical representation for the general solution of the Vekua equation 

3 The two-dimensional Electrical Impedance Equa- 
tion 

Let us consider the two-dimensional case of tlie equation ([T]): 

V • (crVw) = 0. 

As it has been shown in several previous works (see e.g. [5] and [TO])) if f can 
be expressed by means of a separable-variables function 

cr{x,y) = ai{x)a2{y), 



5 



introducing the notations 



P 



yLL, 

0-1 



the equation ([I]) will turn into the Vekua equation 



d^W - —W = 0, 
P 



for which the functions 



Fo = P, Go 



I 

P 



(15) 



(16) 



(17) 



conform a generating pair. From ( 15 1 it is easy to see that this pair is embedded 
into a generating sequence with period k = 2, because p is separable- variables, 
according to Theorem [2j 

3.1 A complete set for boundary value problems of the 
Electrical Impedance Equation 

The possession of an explicit generating sequence, allows the construction of 
the formal powers ( 12 ), so we can approach any solution for ( 16 1, which will be 
closely related with the solutions of ([I]) according to the relations (15). 

As a matter of fact, one special and very important relation between the 
solutions of ([T]) and of (16) was elegantly posed in [5 , and this idea will play a 
central role in the present work. 

Theorem 5 |^ Let us consider the set of formal powers 
|4")(1,2;0), Z(")(^,z;0)^ , 

I J n— 



corresponding to the generating pair Jj7[ ), and let S be a bounded domain 
such that € but ^ F. Then the set of functions valued on T: 



fi?e4")(l,z;0)|r, ReZ^^\i,z-0)\X° , 

) n— 



(18) 



conform a complete system for approaching solutions of the Dirichlet boundary 
value problem corresponding to 

The last statement implies that any boundary condition u\y, imposed for 
the solutions u of ([T]), can be approached according to the expression: 

JV 

lim u\r - V f a'Re^"^ (1, z; 0) |r + a"ReZ^"^ {i, z; 0) |r) = 0, 

n=0 



where a' and a" are real constants. 
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Summarizing, when a separable-variables conductivity function a is given 
within a bounded domain £7, and a boundary condition u|r is imposed for the 
solution o f ([T| ), it will be always possible to construct a finite set of functions, 
subset of (181, such that 



N 



r - XI a'R-e^o"^ (1' z\ 0) |r + a" ReZ^"^ (i, z; 0) |r 



\ri=0 

where e > and I G F. 



dl < e, (19) 



3.2 Construction of a piecewise separable-variables func- 
tion 

One of the main objectives of this work is to contribute into the construction 
of a new theory for the Electrical Impedance Tomography problem. Hence it 
is natural to search for the mathematical tools that will allow us to fully apply 
the modern Pseudoanalytic Function Theory into the analysis, e.g., of medical 
images. 

This means that it is necessary to introduce interpolation methods that, 
given a set of conductivity values defined into a bounded domain on the plane, 
can reach separable-variables functions. One of the first approaches in this 
direction was posed in |10j . and it was properly analysed in jllj. Indeed, the 
last reference is completely dedicated to prove the following assessment. 

Conjecture 1 ///'/ Let a be a function defined within a bounded domain Q £ 
M^, possessing discontinuities only of the first kind. Then it is possible to ap- 
proach a by means of a piecewise separable-variables function of the form: 

■ fniv) ■■ X e [x(M-i),X(M)]. 

(20) 

This separable-variables function can be employed for numerically constructing a 



^ pw \ 



finite set of formal powers of the form { 18), in order to approach solutions for the 



Dirichlet boundary value problem of the two-dimensional Electrical Impedance 
Equation |I]), in the sense of (19), when a boundary condition u\y is imposed. 



A simplified illustration of the procedure for the construction of such piece- 
wise function, is provided along in Figure 1, and a brief explanation of the 
central ideas will be provided now. 

For simplicity, we will consider hereafter the domain VL as the unitary circle, 
noticing that all postulates can be extended to a wide set of bounded domains, 
that are of special interest in many branches of Applied Mathematics, Physics, 
and Engineering. 
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(a) A circular domain sec- (b) Interpolating functions 
tioncd. f{y). 



Figure 1: Simplified illustrations of the steps for approaching piecewise 
separable- variables conductivity functions. 



The first step is to divide the domain into a number M of subsections, 



as shown in the Figure 1(a) Along a strait line crossing every subsection, 
parallel to the y-axis, we will collect a number q of conductivity values, in 
order to introduce an interpolating function / depending only upon the spatial 
variable y. We will posses then a set of M interpolating functions, one for every 
subsection, as shown in Figure l(b)| 



Finally, we state that the conductivity function in every subregion will have 
the separable-variables form 



where j is the number of the subregion, Xj is the common x-coordinate of all 
points collected in the j-subregion, and Kj is a positive real constant such that 
x + Kj ^0, X ^n. 

The postulate remains a conjecture because not any formal extension of the 



theorem posed in 15], about the completeness of the set (18), is known for the 
case of piecewise separable- variables functions within bounded domains. 

Yet, the full set of examples presented in TT] shows that the representation 
( 20 ) is useful for solving the Dirichlet boundary value problem of ([!]) , since the 
numerical calculations succeed to approach the boundary conditions provided by 
known exact solutions, for a variety of conductivities that were both separable- 
variables and non separable-variables by definition. 

This property will be now extended for analysing conductivities upcoming 
from geometrical distributions, precisely as those employed in physical measure- 
ments [in]. To achieve this, we will consider a limit case of the Conjecture [l] 
that will result specially useful for these classes of conductivities. In behalf of 
simplicity, will be considered as the unitary circle, but it is possible to ver- 
ify that the next statements can be generalized for a wide variety of bounded 
domains. 
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Proposition 1 Every conductivity function a, defined within a hounded domain 
J7 G , and possessing only discontinuities of the first kind, can be considered 
the limit case of a piecewise separable-variables conductivity function apw of the 
form (20), when the number M of subsections, and the number q of collected 
values at every subsection, tend to infinity. This is 



a{x,y)= lim (7pw{x,y). (21) 



Moreover, 



X + Ki 



lim ' ^ ^ = 1, 

Af-j-oo Xj + Kj 



where j is the number of the subsection, it follows from Theorem that the 
generating sequence corresponding to this limit case, employed for numerically 
approaching the formal powers, will be periodic with period k = 1. 

Proof 1 Let us consider that Q has been sectioned into M subdomains {f^j}*^]^, 
by employing the set of equidistant y-axis parallel lines 

2j \ 

^3 = ■■ x^f) = -1 + ^; i = 0, 1, An , 

Since coincides with the unitary circle, it can be described by the union of M 
subdomains, defined according to the expressions 

n^+i = {{x,y) : a: e < 1} , j = 0, 1,2, ...,M - 2; 

= {{x,y) : X e [x(M-i),xi^M)], + y'^ <^ ■ 

Let us consider also the set of lines 

^j+i = Xj+1 ■■ Xj+i = X(f) + ^ J = 0, 1, M - 1 



It is clear that 



lim - = 0, j = 0, 1,...,M; 



and in consequence, for^x S Qj', j = 1,2, ...,M; we will have that 

X Xj- 

It immediately follows that 

lim ^±:^ = i. (22) 

M^oo X] + K] 

This implies that every Vtj will be confined into a segment of a parallel line to 
the y-axis. 

Indeed, when the M oo, f2 can be represented by the set of the line seg- 
ments 

{V3 = C : -1 < C < 1, +?;^ < 1} . (23) 
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Beside, since the number q of conductivity values, collected along the line 
ifj+i within every Qj, also tends to infinite, not any interpolation method will 
be required for obtaining the fj{y) functions. They will simply coincide with the 



original values of the function a over the elements of the set (23) 



Finally, from (22), it follows that 



M— >oo \j -\- J\j 



which, according to the last sentence of the Theorem [1| will provoke that the 
corresponding generating sequence, employed to numerically approach some of 
the formal powers will be periodic with period fc = 1. 



4 Numerical solutions for the two-dimensional 
Electrical Impedance Equation 

We now analyse a selected set of analytic conductivity functions for which exact 
solutions are known, in order to examine the effectiveness of the method posed in 
Proposition [T] Thereafter, we will consider conductivity distributions upcoming 
from geometrical distributions, whose analytical representation is, in general, 
unknown, imposing certain boundary conditions that will help us to appreciate 
the behaviour of the technique in this special and important cases. 

A detailed description of the numerical methods used to approach the solu- 
tion, that fulfil the boundary condition of every example, can be found in [3]. 
Since we are considering the unitary circle, and taking into account the valid- 



ity of the expression (19), that as a matter of fact is a Lebesgue integral-type 



operator, let us introduce an inner product for the elements of the finite set 

|Re4"^(l,^;0)|r, ReZ("'(z,z;0)|r}'^ , (24) 
according to the formula 

Re4"'^|r,ReZ^"^V^ =^Re4"^)(0 •ReZ^"^^(/)dZ, 

where I G F, and ni,n2 = 0,1,2,... It follows that we can obtain a set of 
2N + 1 orthonormal functions {uq.}^^q ^, such that we can approach an imposed 
boundary condition it|r according to the expression 

2N+1 

where ba are real constant coefficients. Notice that the apparent lose of one 
base function is because, by virtue of the Definition [5) we have 

Z(°)(*,0;z) = ^, 
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-^1 







■2 

1 



Figure 2: ct = (2 + cos7ra;)(2 + sinTry) within the unitary circle. 



thus IleZQ^\i,0] z) = 0. It should be also noticed that the orthonormalization 
procedure has been performed considering, first, the subset of iV + 1 functions 

|Re4")(l,0;z)r , (25) 

I. J n— 

followed by the subset of functions 

|ReZ^")(z,0;z)}'^ . (26) 

L J n— 1 

This remark is important for adequately examining the illustrations where the 
absolute values of the coefficients are displayed. 

4.1 The sinusoidal case. 
Proposition 2 Let 

a = (2 + cos lux) (2 + sin ujy) . (27) 

Then the function 

2 /tan^^\ 2 /l + 2tan5f\ ^ ^ 

. = ^ arctan (^-^ j + ^ arctan ^ ) ; (28) 

is a particular solution of the Electrical Impedance Equation 

This case was selected because of the variations of conductivity that take place 
within the unitary circle. Unfortunately, for avoiding the indetermination of 
the tangent functions contained into the particular solution ( 28 1 , we could only 
consider the case a; — > tt. An illustration of this conductivity is given in the 
Figure 2. 

The number of points located around F, when imposing the boundary con- 
dition, does not necessarily has to coincide with the number a of orthonormal 
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Table 1: Values of the coefficients corresponding to the boundary value 
problem with cr = (2 + C0S7rx)(2 + sin Try). 



bo 


bi 


62 


64 


bi9 


620 


622 


623 


9.022 


30.180 


12.446 


3.914 


-38.846 


4.632 


2.820 


1.627 



functions, as it is shown in [6^, where the collocation method is employed. Nev- 
ertheless, on behalf of simplicity, this work will consider 35 equally distributed 
points around the perimeter of the unitary circle, and a = 35 base functions Ua- 
The boundary condition will be obtained by evaluating the solution m, presented 



in (28 1, over this set of points. 



The absolute error £ is defined as the classical Lebesguc norm 

£ = U Lil) - bo^ua^m dl] ; (29) 



where u{l) represents the solution (28) valued on the boundary F, and N = 17. 
The result of this integral will be approached using the standard trapezoidal 
method over 1000 equally distributed points on the segment [0, 27r]. 



The Figure 3(a) displays a logarithmic plot of the absolute magnitudes cor- 
responding to every coefficient ba used for approaching the boundary condition. 
We shall remember that the first 18 coefficients correspond to orthonormal sys- 
tem obtained from the set whereas the remaining 17 correspond to the set 
([26|. The Figure [3(b)] displays a comparison between the imposed boundary 



condition ( 28 ) valued on F, drew in blue, and the approached solution, in red. 
It is not possible to detect any difference at first sight in this case. Beside, the 
obtained error £ — 0.155 x lO^'^, thus it is possible to assess the approach is 
adequate. 

The Table [l] show the numerical values of some of the most relevant coef- 
ficients ba used in the reconstruction. It is important to pay attention to the 
number a selected, since, in general, they are not presented in consecutive order. 

4.2 The Lorentzian cases. 

Proposition 3 111}/ Let the conductivity function be 

' > ^ (30) 



(x-isy + 0.1 1 Vr + 0.1 



Thus the function 

3 

is a solution of 



+ 0.1(x-/3 + y), (31) 



12 




13 




(c) Case when /3 = 1. 

Figure 4: Illustrations for the case a = {{x - jS^ + 0.1)"^(y^ + 0.1)~^ 
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Table 2: Values of the coefRcients corresponding to the boundary value 
problem with a = {x'^ + 0.1)"^ (?/^ + 0.1)"^ 



bi 


63 


^'5 


bg 


big 


621 


^^23 


b25 


7.923 


0.563 


-1.126 


-0.261 


-7.725 


2.220 


-0.331 


0.058 



Table 3: Values of the coefficients ba corresponding to the boundary value 
problem with a = {{x - 0.5)^ -|- 0.1)-'^ {y^ + 0.1)-^. 



bo 


61 


b2 


^3 


619 


&20 


&21 


&22 


3.743 


15.051 


6.542 


-3.125 


-7.723 


2.160 


0.578 


-0.179 



/3 = 0, 0.5, 1. The Figures [4(a)l [4(b)] and \A{c] 



Three cases will be considered: 
display each example. 

Following the same logical steps of the previous subsection, the Figures 5(a) 
|5(b)| and 5(c)[ show logarithmic plots of the absolute values of the coefficients 
ba corresponding to the cases /? = 0, 0.5, 1; respectively. The illustrations with 
the comparisons between the boundary conditions obtained from (31) and the 
approached solutions were omitted, because once more, it was not possible to 
notice any difference between the pairs of curves. 

The absolute error for the case /3 = was £ = 0.592 x 10~'^, the one cor- 
responding to /? = 0.5 was £ — 1.4 x 10"'^, and such belonging to /3 = 1 was 
f = 2.9 X 10-3. 

The Tables [2] [3] and [4] show some of the most relevant coefRcients corre- 
sponding to the approached solutions of the cases (3—1, 0.5, 1. 



4.3 Conductivities corresponding to geometrical distribu- 
tions. 

Let us consider first a piecewise conductivity function, in polar coordinates, of 
the form: 



100 : 


r G [0,0.2); 


30 : 


r e [0.2,0.4); 


20 : 


r e [0.4,0.6); 


15 : 


r e [0.6,0.8); 


30 : 


r G [0.8,1). 



(32) 



Table 4: Values of the coefficients ba corresponding to the boundary value 
problem with a = {{x - 1)^ + 0.l)-'^{y'^ + 0.1)"^ 



bo 


bi 


b2 


^3 


&4 


&5 


be 


619 


6.355 


37.838 


11.3959 


-10.118 


-3.517 


-3.311 


-1.815 


-7.856 



15 
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(a) Absolute values of the 35 coefficients employed for approaching the boundary 
condition when 13 = 0. 



10' 
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(b) Absolute values of the 35 coefficients employed for approaching the boundary 
condition when f} = 0.5. 



(c) Absolute values of the 35 coefficients employed for approaching the boundary 
condition when 13 = 1. 
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Figure 5: Coefficients for the cases cr = ((a; - /S)^ + 0.iy^{y'^ + 0.1)"^; /3 
1, 0.5, 1. 
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Figure 6: Conductivity defined according to the expression (32 1 
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Figure 7: Conductivity defined according to the expression (32) 



Here r denotes the radio. The Figure |??| iUustrates this conductivity. To select 
a boundary condition for the geometrical cases, without performing physical 
measurements, it is not a trivial task. Nevertheless, the conductivity defined in 



(32) shall be somehow related with the Lorentzian cases previously studied. For 



this reason, the boundary condition will be precisely the expression ( 32 ) 



+ 0.l{x + y). 



The Figure [??] displays the absolute value of the coefficients 6^. Surprisingly, 
only four coefficients were significant for fulfilling the imposed boundary condi- 
tion, and the total error was £ = 1.486 x 10"^''. Again, the graphic comparing 
the boundary condition and the approach is omitted, since not any difference is 
visible. 

The Table [S] presents the values of the four most significant coefficients for 
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Table 5: Values of the coefficients 6^ corresponding to the boundary value 



problem with cr defined in ( 32 ) 



bi 


^'3 


bi9 


b2i 


7.826 


1.863 


-7.826 


1.863 



this experiment, which also exhibit a very interesting symmetry. 

Let us consider now the conductivity functions illustrated in the Figures 
a)] |8(b)| and 8(c) The blue sections represent conductivity values of a 



10, whereas the red circles posses a = 100. The red disk in the Figure 8(a) 
corresponds to the equation 

< 0.2, 

the red disk in the Figure |8(b)| is traced according to 

(x- 0.6)2 < Q 2, 



whereas the red disk in the Figure 8(c) is given by 



{x~0.79f < 0.2. 



The boundary condition for the conductivity illustrated in the Figure |8(a) 
will be 

u=—^+0.1{x + y), 



the condition for the conductivity function showed in the Figure [8(b) [ will be 
u= ^^~'f + y\ 0.lix^0.6 + y); 



and finally, the condition imposed for the conductivity of the Figure 
selected as 



was 



(x~ 0.79)3 + r 



+ 0.1(2; -0.79 + 2/). 



The Figures 9(a)[ |9(b)| and 9(c) display the logarithmic graphics of the 
absolute values belonging to the coefficients ba- Once more, it is interesting 
that the quantity of significant values is relatively small. The absolute error 

1.623 X 10-1-* for the case 



approached for every case were: £ 



£ 



7.3 X 10-3 for the case |8(b) and f = 3.9 x lO'^ for the case |8(c) 

As before, the Tables [6j [7] and [8] contain the values of the most significant 
coefficients for each of the last three examples. 



4.4 A triangular surface within the unitary circle. 

This is the last and one of the most interesting cases analysed in this work. We 
consider a triangle within the bounded domain, as displayed in the Figure 10, 
which somehow could be considered an interesting challenge when solving the 
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(a) Disk ^2 + < Q 2 



(b) Disk {x - 0.6)2 + < 0.2 




-1 1 

(c) Disk (x - 0.79)2 + y'^ < 0.2 



Figure 8: Conductivity functions with a single disk within the domain. The 
red surfaces represent a — 100, whereas the blue sections denote a — 10. 



Table 6: Values of the coefficients ba corresponding to the boundary value 
problem with a posed in 8(. 



bi 


&3 


bi9 


b2i 


7.826 


1.863 


-7.826 


1.863 



Table 7: Values of the coefficients ba corresponding to the boundary value 



problem with a posed in 8(b) 



^^0 


bi 


b2 


^^3 


^19 


b2i 


-13.655 


15.883 


-6.712 


1.830 


-7.838 


1.819 



Table 8: Values of the coefficients ba corresponding to the boundary value 



problem with a posed in 8(b) 



bo 


bi 


b2 


b3 


^19 


b2i 


-19.585 


21.781 


-8.832 


1.861 


-7.830 


1.846 
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(a) Absolute values of the 35 coefficients employed for approaching the boundary 
condition of the conductivity |8(a)[ 



10' 



(b) Absolute values of the 35 coefficients employed for approaching the boundary 
condition of the conductivity |8(b)] 



(c) Absolute values of the 35 coefficients employed for approaching the boundary 
condition of the conductivity |8(c)| 
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Figure 9: CoefRcients ba for the cases of disks within the unitary circle. 
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Figure 10: Conductivity containing a triangular figure. The red area represents 
a — 100, whereas the blue is cr = 10. 

Dirichlet problem of the two-dimensional Electrical Impedance Equation ([T]), 
using numerical methods based upon variations of the Finite Element Method, 
one of the finest known numerical tools for this kind of problems. 

The presence of corners requires special attention. Because of this, at least 
one radial trajectory in the calculations of the formal powers was force to cross 
over every corner of the triangle. Yet, the selection of the boundary condition is 
a completely different question. Indeed, this case will barely allow the imposition 
of theoretical conditions reaching acceptable convergence results. 

For this reason, without denying the arbitrary selection of the function, the 
boundary condition was established as 



This is indeed the only case for which 32 formal powers were numerically ap- 
proached, reaching a system of 61 orthonormal functions. The total error was 
£ = 0.057, but the graphic containing the comparison between the boundary 
condition and the approached solution, traced in Figure 11, deserves special 
attention. 

Let I G [0, 27r). The highest divergence is located around the values I = \ 
and Z = as it is pointed out in the corresponding decimal representation 
on the Figure 11. These coincide precisely with the radii where the corners of 
the triangle are located, which posses interesting questions. Nevertheless, it is 
convenient to remember that the boundary condition was arbitrary selected, 
thus perhaps it is better to postpone the discussions until the moment that 
physical measurements provide the boundary conditions to be approached. 

As for the rest of the posed examples, the Table |9] contains some of the most 
representative values of the coefficients used for approaching the solution. 



{x - 0.6)3 ^ y3 
3 



+ 0.1 (x- 0.6 + 2/). 
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-1 

-1.2 
-1.4 











\ — Boundary condition 
\ —Approached solution 

















































0.786 2 4 5.492 6 



Figure 11: Comparison between the boundary condition and the approached 
solution of the conductivity function displayed in the Figure 10. 

Table 9: Values of the coefficients ba corresponding to the boundary value 
problem with a posed in Figure 10. 



bo 


bi 


&2 


bs 


^32 


^21 


-13.617 


15.919 


-6.281 


2.929 


-7.845 


1.231 



5 Conclusions 

On of the most important contributions of the present work is the full opening 
of the path for applying the Modern Theory of Pseudoanalytic Functions, into 
the analysis of wide class of conductivity functions upcoming from physical 
problems, by virtue of the Proposition [T] 

Tested in a variety of examples, the numerical analysis based upon the cited 
proposition, succeed to approach the imposed boundary conditions with con- 
siderable accuracy, presenting the highest divergence only in a kind of problem 
that is well known for its complexity when analysing boundary value problems 
for the two-dimensional Electrical Impedance Equation. 

An immediate implication of these results is the possibility of analysing 
most classes of images corresponding to classical applications of the Electrical 
Impedance Tomography, as it is the Medical Imaging clinical monitoring. In this 
precise direction, it is possible to assure that the search for patterns of change 
in the boundary electric potentials, when changes in the conductivity within 
the domain of interest are taking place, is completely viable by applying the 
techniques posed before. It is important to remark that all presented results can 
be extended without mayor complications to a wide class of bounded domains, 
beside the unitary circle, that on behalf of simplicity was considered in these 
pages. 

Indeed, most problems of Mathematical Physics, closely related with the 
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Modern Pseudoanalytic Function Theory, as those elegantly studied in [5], as 
well as in |7j, where a special case of the Fokker-Planck equation is posed, could 
well be susceptible for this kind of analysis. 

Still, it is the Electrical Impedance Tomography problem one of the most 
interesting fields for applying the techniques suggested in this work. 
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